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Abstract 

A kinematic model based on the superposition of p + p collisions, relativistic geometry and 
final-state hadronic rescattering is used to predict various hadronic observables in y / Siv r iv = 2.76 
TeV Pb + Pb collisions. Predictions for rapidity and transverse momentum distributions, elliptic 
flow, and two-boson femtoscopy are presented. A short proper time for hadronization is assumed. 
Previous calculations using this model which were performed for y / sjv r /v = 200 GeV Au + Au 
collisions were shown to describe reasonably well the trends of observables in experiments carried 
out at that energy, giving the present predictions for Pb + Pb at higher energy some degree of 
credibility. 
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I. INTRODUCTION 



The present plan for running at the CERN Large Hadron Collider (LHC) calls for Pb+Pb 
collisions at ^snn = 2.76 TeV to commence in November, 2010. These will be the highest 
energy heavy ion collisions ever produced in the laboratory, and it is hoped that exotic 
phenomena will be observed which will expand our knowledge of the properties of highly 
excited matter jl]. 

In this context, it is the goal of the present paper to make predictions for common hadronic 
observables which will be initially measured by LHC experiments in these high-energy Pb + 
Pb collisions. Predictions for spectra (i.e. rapidity and transverse momentum distributions), 
elliptic flow, and two-boson femtoscopy are presented. A hadronic rescattering model in 
which the initial state is determined by the superposition of proton-proton collisions has 
been chosen to make these predictions. The advantages of this model for this purpose are 
1) the model has been shown to describe the overall trends of hadronic observables in lower 
energy Au + Au collisions at ^snn = 0.20 TeV from the Relativistic Heavy Ion Collider 
(RHIC) 2], and 2) the model is easily scalable to LHC energies. These will be "limiting 
case scenario" predictions in the sense that only hadrons are used as the degrees of freedom 
in this model even at the early stages of the collision where parton degrees of freedom are 
thought to be more appropriate, i.e. a short proper time for hadronization is assumed. 

The paper is organized into the following sections: Section [III gives a brief description of 
the model, Section III II presents predictions from the model for ^s^n — 2.76 TeV Pb + Pb 
collisions, and Section IIVI gives a summary and conclusions. 



II. DESCRIPTION OF THE MODEL 

The model calculations are carried out in five main steps: A) generate hadrons in p + p 
collisions from PYTHIA, B) superpose p+p collisions in the geometry of the colliding nuclei, 
C) employ a simple space-time geometry picture for the hadronization of the PYTHIA- 
generated hadrons, D) calculate the effects of final-state rescattering among the hadrons, 
and E) calculate the hadronic observables. These steps will now be discussed in more detail. 
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A. Generation of the p + p collisions with PYTHIA 

The p + p collisions were modeled with the PYTHIA code [3], version 6.409. The internal 
parton distribution functions "CTEQ 5L" (leading order) were used in these calculations. 
Events were generated in "minimum bias" mode, i.e. setting the low-pr cutoff for parton- 
parton collisions to zero (or in terms of the actual PYTHIA parameter, ckin(3) = 0) and 
excluding elastic and diffractive collisions (PYTHIA parameter msel = 1). Runs were made 
at y/s = 2.76 TeV to simulate the upcoming LHC collisions. Information saved from a 
PYTHIA run for use in the next step of the procedure were the momenta and identities of 
the "direct" (i.e. redundancies removed) hadrons (all charge states) ir, K, p, n, A, A, p, 
u, rj, rj', <f), and K*. These particles were chosen since they are the most common hadrons 
produced and thus should have the greatest effect on the hadronic observables in these 
calculations. 

B. Superposition of p + p events to simulate heavy-ion collisions 

An assumption of the model is that an adequate job of describing the heavy-ion collision 
can be obtained by superposing PYTHIA-generated p + p collisions calculated at the beam 
yfs within the collision geometry of the colliding nuclei. Specifically, for a collision of impact 
parameter b, if f(b) is the fraction of the overlap volume of the participating parts of the 
nuclei such that f(b = 0) = 1 and f(b = 2R) = 0, where R = 1.2A 1 / 3 and A is the mass 
number of the nuclei, then the number of p + p collisions to be superposed will be f{b)A. 
The positions of the superposed p + p pairs are randomly distributed in the overlap volume 
and then projected onto the x — y plane which is transverse to the beam axis defined in the 
^-direction. The coordinates for a particular p + p pair are defined as x pp , y pp , and z pp = 0. 
The positions of the hadrons produced in one of these p+p collisions are defined with respect 
to the position so obtained of the superposed p + p collision (see below). 

As was done in similar calculations for RHIC collisions to give better agreement with 
experimental dn/drj distributions^, a lower multiplicity cut was applied to the p+p collisions 
used in the present calculations which rejected the lowest 20% of the events. The spirit 
of this cut is to partially compensate for the fact that there is no reinteraction of primary 
nucleons from the projectile-target system in this model. 
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C. The space-time geometry picture for hadronization 

The space-time geometry picture for hadronization from a superposed p + p collision 
located at (x pp , y pp ) consists of the emission of a PYTHIA particle from a thin uniform 
disk of radius 1 fm in the x — y plane followed by its hadronization which occurs in the 
proper time of the particle, r. The space-time coordinates at hadronization in the lab frame 
(xh,Vh, Zh,th) for a particle with momentum coordinates (p x ,P y ,Pz), energy E, rest mass 
mo, and transverse disk coordinates (xo,yo), which are chosen randomly on the disk, can 
then be written as 



x h = x pp + x + t— (1) 
m 

D 

Vh = y PP + yo + r— (2) 

Zh = T— (3) 

m 

th = r— (4) 
m 

Eqs. (1) and (2) show the initial expansion in the transverse direction now present in 
the model. The simplicity of this geometric picture is now clear: it is just an expression 
of causality with the assumption that all particles hadronize with the same proper time, r. 
A similar hadronization picture (with an initial point source) has been applied to e + — e~ 
collisions [4]. For all results presented in this work, r will be set to 0.1 fm/c as was done in 
applying the present model to calculating predictions for RHIC Au + Au collisions [2] and 
Tevatron p + p collisions pf • 



D. Final-state hadronic rescattering 

The hadronic rescattering calculational method used is similar to that employed in pre- 



vious studies 



BE 



7j . Rescattering is simulated with a semi-classical Monte Carlo calculation 
which assumes strong binary collisions between hadrons. Relativistic kinematics is used 
throughout. The hadrons considered in the calculation are the most common ones: pions, 
kaons, nucleons and lambdas (it, K, N, and A), and the p, u, rj, rf, <ft, A, and K* resonances. 
For simplicity, the calculation is isospin averaged (e.g. no distinction is made among a 7r + , 
7t°, and 7r~). 



The rescattering calculation finishes with the freeze out and decay of all particles. Starting 
from the initial stage (t — fm/c), the positions of all particles in each event are allowed to 
evolve in time in small time steps (At = 0.5 fm/c) according to their initial momenta. At 
each time step each particle is checked to see a) if it has hadronized (t > t^, where th is given 
in Eq. (4)), b) if it decays, and c) if it is sufficiently close to another particle to scatter with 
it. Isosp in- averaged s-wave and p-wave cross sections for meson scattering are obtained from 
Prakash et al. jsj and other cross sections are estimated from fits to hadron scattering data 
in the Review of Particle Physics Both elastic and inelastic collisions are included. The 
calculation is carried out to 400 fm/c which allows enough time for the rescattering to finish 
(as a test, calculations were also carried out for longer times with no changes in the results). 
Note that when this cutoff time is reached, all un-decayed resonances are allowed to decay 
with their natural lifetimes and their projected decay positions and times are recorded. 

The rescattering calculation is described in more detail elsewhere js, 7]. The validity of 
the numerical methods used in the rescattering code have been studied using the subdivision 
method, the results of which have verified that the methods used are valid [10]. 



E. Calculation of the hadronic observables 

Model runs are made to be "minimum bias" by having the impact parameters of collisions 
follow the distribution da/db oc b, where < b < 2R. Observables are then calculated from 
the model in the appropriate centrality bin by making multiplicity cuts as normally done 
in experiments, as well as kinematic cuts on rapidity and pr- For the present study, a full- 
calculation 3200 event minimum bias run was made from the model for y / sjv r /v = 2.76 TeV 
Pb + Pb collisions which was then used to calculate all of the hadronic observables shown. 
In addition, a 3200 event minimum bias run with rescattering turned off in the model was 
also made for comparison to study the importance of rescattering in these observables. 

III. PREDICTIONS FOR ^/IJ^ = 2.76 TEV Pb + Pb COLLISIONS 

As mentioned earlier, predictions for the hadronic observables spectra, elliptic flow, and 
two-boson femtoscopy have been made with the present model for ^/snn = 2.76 TeV Pb+Pb 
collisions. Results for each of these observables are presented separately below. 
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FIG. 1. Model charged-hadron pseudorapidity distributions for Pb + Pb collisions at ^/snn = 2.76 
TeV for minimum bias and — 5% centrality bins. Also shown is a comparison of non-single 
dim-active PYTHIA p + p collisions at yfs = 2.76 TeV with a measurement from ALICE scaled to 
this energy. 



A. Spectra 

Figures [U, [2} and [3] show predictions from the present model for charged hadron spectra. 

Figure Q] shows charged-hadron pseudorapidity distributions for Pb + Pb collisions at 
t/snn = 2.76 TeV for minimum bias and — 5% centrality bins. Also shown is a comparison 
of non- single diffractive PYTHIA p + p collisions at i/i = 2.76 TeV with a measurement 
of non-single-diffractive dn/drj from the LHC ALICE experiment made at \fs = 2.36 TeV 



11] , approximately scaled to the slightly higher energy by y^. 76/2. 36. Since the model is 
"isospin averaged", the model distributions are multiplied by 2/3 to approximate all charged 
particles. Looking at the p+p distribution first, it is seen that PYTHIA is in good agreement 
with the scaled ALICE measurement at mid-rapidity, giving additional confidence in using 
PYTHIA p + p collisions at this energy in the superposition for the Pb + Pb collisions. 
Looking at the distributions for Pb + Pb, the mid-rapidity dn/drj values for minimum bias 




FIG. 2. Model charged-hadron pt distributions for the same conditions as shown in Figure [TJ 



and central collisions are found to be 323 and 1174, respectively. For central collisions, this 



is about twice the value found in RHIC ^s^n = 200 GeV Au + Au collisions [12| , and it is 
at the lower end of the range of predictions which have been recently made of 1500 — 4000 in 
central collisions using various extrapolations of RHIC experimental rapidity densities []] . 

Figure |2] shows charged-hadron distributions for Pb+Pb collisions at a/saT/v — 2.76 TeV 
for minimum bias and 0—5% centrality bins. Also shown for comparison is the distribution 
for non-single diffractive PYTHIA p + p collisions at yfs = 2.76 TeV. To approximate 
(h + + h~)/2 for the model, the model distributions are multiplied by 1/3. The — 5% 
centrality charged particle px distribution for pt > 5 GeV/c is predicted to be about two 
orders of magnitude larger at the LHC compared with RHIC [13]. This is an expected 
consequence of the higher y^sJvjv in the LHC collisions that the px distributions at high p T 
should be greatly enhanced. 

Figure [3] shows the ratio of minimum bias Pb + Pb to p + p pt distributions from Figure 
[2] (indicated as "Rescattering ON"). Also shown for comparison is the similar ratio with 
rescattering turned off in the model. This ratio is essentially the "unnormalized" Raa 
observable 14j. Studying the high p T behavior of Raa is thought to be a way of more directly 
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FIG. 3. Ratio of minimum bias Pb + Pb to p + p pt distributions from Figure [2] (indicated as 
"Rescattering ON"). Also shown is the similar ratio with rescattering turned off in the model. 



studying QCD processes, such as jets, in heavy-ion collisions. Since the present model is 
based on using PYTHIA which uses QCD processes in calculating p + p collisions, the model 
should contain these effects and thus should be suitable for comparing with experiments 
which measure these observables. As seen in Figure [3j for the full rescattering calculation 
the ratio is suppressed for pt > 2 GeV/c compared with the case of rescattering turned off. 
This high-p^ suppression is similar to that observed in RHIC collisions jl^. In the present 
model calculations this suppression is clearly seen to be due to hadronic rescattering. 
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B. Elliptic flow 



The elliptic flow variable, V2, is defined as 

V 2 =< cos(20) > (5) 
4> = arctanf— ) 

Px 

where "<>" implies a sum over particles in an event and a sum over events and where 
p x and p y are the x and y components of the particle momentum, and x is in the impact 
parameter direction, i.e. reaction plane direction, and y is in the direction perpendicular to 
the reaction plane. The V2 variable is calculated from the model using Eq. (JSJ) and taking 
the reaction plane to be the model x — z plane. 

Figure H] shows model V2 vs. px plots for t/snn = 2.76 TeV Pb + Pb collisions for 
minimum bias and — 5% centrality bins, all hadrons, and — 1 < 77 < 1. Also shown is 
the minimum bias case with rescattering turned off. For the minimum bias case with the 
full calculation, V2 is seen to increase with increasing px, peaking at a value of about 0.14 
and then decreasing for px increasing beyond 2.5 GeV/c. For the — 5% centrality case, 
V2 is also seen to increase with increasing px and to have a significantly smaller magnitude 
than the minimum bias case. These general behaviors of V2 are also observed in RHIC 



collisions 
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16]. When rescattering is turned off, V2 — > for all px, demonstrating that 



hadronic rescattering accounts completely for the elliptic flow signal in this model. 



C. Two-boson femtoscopy (Hanbury-Brown-Twiss interferometry) 

Figures |5||9] and Tables Ul lIIII show predictions from the model for two-pion and two-kaon 

n 

HBT for t/snn = 2.76 TeV Pb + Pb collisions. For the HBT[17] calculations from the model, 
the three-dimensional two-boson correlation function is formed and a Gaussian function in 
momentum difference variables is fitted to it to extract the boson source parameters. Boson 
statistics are introduced after the rescattering has finished (i.e. when all particles have 
"frozen out" ) using the standard method of pair- wise symmetrization of bosons in a plane- 
wave approximation [18]. The three-dimensional correlation function, C \Q S ide,Q out, Q 'long), 



is then calculated in terms of the momentum-difference variables Q S ide, which points in the 
direction of the sum of the two boson momenta in the transverse plane, Q ou t, which points 
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FIG. 4. Model V2 vs. pr plots for + Pb collisions for minimum bias and — 5% centrality bins, 
all hadrons, and — 1 < 77 < 1. Also shown is the minimum bias case with rescattering turned off. 



perpendicular to Q S ide in the transverse plane and the longitudinal variable along the beam 
direction Qi ong . 

The final step in the calculation is extracting fit parameters by fitting a Gaussian param- 
eterization to the model-generated two-boson correlation function given by, 

C^Qsidey Qouty Qlong) 

G[l + Aexp(— Q side Rl ide — Ql ut Rlut ~ QlongRlong)) (6) 

where the i?-parameters, called the radius parameters, are associated with each momentum- 
difference variable direction, G is a normalization constant, and A is the usual empirical 
parameter added to help in the fitting of Eq. ([6]) to the actual correlation function (A = 1 

ongitudinally 



in the ideal case). The fit is carried out in the conventional LCMS frame 
comoving system) in which the longitudinal boson pair momentum vanishes [19|]. Figure [5] 
shows a sample projected two-pion correlation function from the model with projected fit 
to Eq. ([6]). It is seen that Eq. ([6]) fits the model-generated correlation function quite well. 
Figure [6] shows the dependence of the model pion source parameters for Pb + Pb collisions 
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FIG. 5. Sample two-pion correlation function with Gaussian fit projected onto the Q ou t, Qsidei an d 
Qiong axes from the model. The collision centrality is — 5% and the kx bin is 0.1 — 0.2 GeV/c. 



on k~T for centrality — 5% and — 1 < rj < 1 for the full calculation compared with the 
calculation with rescattering turned off. For the full calculation, the radius parameters 
are all predicted to decrease with increasing showing; the effects of "flow" as has been 



19L |20|. The overall scales of R S ide 
IC Au + Au collisions, 



201 ] . The A-parameter 



observed in RHIC Au + Au collisions and elsewhere 
and Rout predicted for Pb + Pb are comparable to those seen in RH] 
whereas Ri ong is predicted to be about 25% larger than at RHIC 
is seen to be mostly independent of with a value of about 0.6, which is significantly 
less than the "ideal HBT case" of A = 1. The main effect causing A < 1 in the model is 
the presence of long-lived resonances such as rj and rj' which decay into pions late in the 
collision thus suppressing the correlation function. Looking at the radius parameters with 
rescattering turned off in Figure El it is seen that their dependence on k T mostly disappears 
and their scales are significantly reduced compared with the full calculation, showing the 
strong influence that rescattering has on the HBT parameters in this model. Comparing 
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FIG. 6. Model pion source parameters vs. kx for Pb + Pb collisions, — 5% centrality, — 1 < r) < 1, 
and with and without rescattering turned on. 



A-parameters, in the case with rescattering turned off the A-parameter is seen to increase 
with increasing kr unlike for the full calculation. An explanation for this behavior is that 
with rescattering turned off, the pion source is less Gaussian than with rescattering turned 
on and thus the A-parameter tries to compensate in the fit for this effect. 

Figure [7] compares the model pion source parameters vs. kx for Pb + Pb collisions, where 
— 1 < f] < 1, for the two centralities — 5% and 62 — 72%. Tables HI and HT1 give the values 
for the plots in Figure [71 The more peripheral centrality case, i.e. 62 — 72%, is seen to have 
some qualitative similarities with the "Rescattering OFF" calculation shown in Figure El in 
that the radius parameters are seen to have weaker dependences on kx than for the more 
central case, i.e. — 5%, and are significantly smaller in magnitude. The A-parameter kx 
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FIG. 7. Model pion source parameters vs. kx for P6 + Pb collisions, — 1 < rj < 1, — 5% and 
62 — 72% centralities. 



dependence for the peripheral case also resembles that for the "Rescattering OFF" case in 
that A more or less increases with increasing as opposed to being mostly independent of 
kx as for the full-calculation central case. The explanation for the similarities between the 
peripheral case and the "Rescattering OFF" case is due to the smaller particle multiplicity 
for the peripheral collisions and thus less rescattering present than for the central case, 
determined by the smaller initial geometric overlap of the projectile-target system in the 
peripheral case. 

Model predictions have also been made for two-kaon HBT as shown in Figures |S] and El 
Figure [8] shows a sample two-kaon correlation function with Gaussian fit projected onto the 
Qout, Qside, and Qiong axes from the model for the collision centrality — 5%, — 1 < r\ < 1 
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TABLE I. Model pion source parameters vs. fcy for y / J/v7V = 2.76 TeV Pb ± Pb collisions, — 1 < 
77 < 1, — 5% centrality. 



Jct range 


< kx > 


A 


R d 


R t 




(GeV/c) 


(GeV/c) 




(fm) 


(fm) 


(fm) 


0.10-0.20 


0.15 


0.613 ± 0.000 


5.33 ±0.00 


6.40 ± 0.00 


8.25 ±0.00 


0.20 - 0.30 


0.25 


0.608 ± 0.007 


4.68 ±0.05 


6.21 ±0.07 


6.65 ±0.06 


0.30 - 0.40 


0.35 


0.592 ± 0.003 


4.40 ±0.07 


5.93 ±0.11 


5.58 ±0.05 


0.40 - 0.50 


0.45 


0.582 ± 0.022 


4.20 ±0.06 


5.74 ±0.21 


4.46 ± 0.06 


0.50-0.70 


0.57 


0.622 ± 0.062 


3.99 ±0.13 


5.10 ±0.29 


4.05 ±0.21 


0.70-1.00 


0.76 


0.552 ± 0.038 


3.63 ±0.23 


3.67 ±0.23 


2.53 ±0.32 


TABLE II. 


Model pion source 


parameters vs. 


k T for yJs NN = 


2.76 TeV Pb + Pb collisions, 


-1 < 77 < 1, 


62 - 72% centrality. 








kx range 


< k T > 


A 


Rside 


Rout 


Rlong 


(GeV/c) 


(GeV/c) 




(fm) 


(fm) 


(fm) 


0.10-0.20 


0.15 


0.364 ± 0.003 


2.20 ±0.02 


2.87 ±0.01 


3.31 ±0.02 


0.20 - 0.40 


0.29 


0.430 ± 0.004 


2.33 ±0.01 


2.91 ±0.02 


2.79 ± 0.02 


0.40 - 0.60 


0.48 


0.647 ± 0.009 


2.24 ±0.06 


2.89 ± 0.05 


2.37 ±0.07 


0.60-1.00 


0.68 


0.564 ± 0.006 


0.88 ±0.03 


2.07 ± 0.03 


1.65 ±0.03 



and fcr bin 0.1 — 1.0 GeV/c. The Gaussian fit is seen to provide a reasonable fit to the 
model correlation function. Figure [9] presents a comparison of model pion and kaon source 
parameters vs. for Pb + Pb collisions, —1 < r\ < 1, centralities — 5% and 62 — 72%. 
Table II III gives the values for the two-kaon source parameters plotted in Figure [9j A large 
^ bin of 0.10 < kx < 1.00 GeV/c with an average hp of < kx >= 0.39 GeV/c was used 
for the two-kaon calculations in order to obtain reasonable statistical errors from the 3200 
event minimum bias Pb + Pb run used in the present study. Also, while a pseudorapidity 
range of —1 < 77 < 1 was used for the — 5% centrality two-kaon calculation, a range of 
—4 < 77 < 4 was used for the 62 — 72% centrality two-kaon calculation in order to obtain 
reasonable statistical errors. As seen in Figure [9], the two-kaon calculations obey "mr 
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FIG. 8. Sample two-kaon correlation function with Gaussian fit projected onto the Q ou t, Qside, and 
Qiong axes from the model. The collision centrality is — 5% and the kr bin is 0.1 — 1.0 GeV/c. 



TABLE III. Model kaon source parameters vs. centrality bin for ^/snn = 2.76 TeV Pb + Pb 
collisions, 0.10 < k T < 1.00 GeV/c (< k T >= 0.39 GeV/c). 



Centrality bin 


A 




Pout 


^long 






(fm) 


(fm) 


(fm) 


0-5% 


0.781 ± 0.030 


3.60 ± 0.24 


4.47 ±0.16 


4.51 ±0.35 


62 - 72% 


0.547 ± 0.054 


2.04 ± 0.26 


2.50 ±0.13 


2.49 ± 0.49 



scaling" reasonably well with the two-pion calculations considering the large kx bin that 
was necessary to use for the kaons. 
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FIG. 9. Comparison of model pion and kaon source parameters vs. my for Pb + Pb collisions, 
-1< T) < 1, - 5% and 62 - 72% centralities. 



IV. SUMMARY AND CONCLUSIONS 



A kinematic model based on the superposition of PYTHIA-generated p + p collisions, 
relativistic geometry and final-state hadronic rescattering has been used in the present work 
to predict several hadronic observables in y / J/v r /v = 2.76 TeV Pb + Pb collisions. A short 
proper time for hadronization of r = 0.1 fm/c has been assumed as in previous studies 
with this model which have shown qualitative agreement with experiments. Predictions for 
rapidity and transverse momentum distributions, elliptic flow, and two-boson femtoscopy 
have been presented which will likely be among the first observables to be extracted from 
analyses of the first Pb + Pb data from the LHC. 
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The most noticeable features of the predictions from the present model study which have 
been presented are summarized below: 

• dn/dr] near mid-rapidity for charged particles in LHC Pb + Pb collisions is predicted 
to be about 1200 for a — 5% centrality window. This puts its value at the lower end 
of the range of predictions which have been recently made of 1500 — 4000 in central 
collisions using various extrapolations of RHIC experimental rapidity densities. 

• The — 5% centrality charged particle px distribution for p T > 5 GeV/c is predicted 
to be about two orders of magnitude larger at the LHC compared with RHIC. 

• High-p T suppression of the Raa is predicted to be present in LHC Pb + Pb collisions 
as has been seen in RHIC collisions. 

• Elliptic flow for charged hadrons is predicted to be comparable to that seen in RHIC 
collisions. 

• Two-pion HBT radius parameters from LHC Pb + Pb are predicted to be comparable 
in scale to those from RHIC Au + Au collisions for R si d e and Rout, and 25% larger for 
Riong and are predicted to show decreasing magnitude with increasing kx, i.e. "flow" 
effects. Also, two-kaon HBT radius parameters are predicted to show "m T -scaling." 

As mentioned earlier, these are "limiting case scenario" predictions from a model that 
only considers hadronic degrees of freedom, i.e. which assumes that only "ordinary" physics 
processes are taking place in these collisions. Of course the "best case scenario" will be that 
these predictions disagree wildly with the actual measured hadronic observables which will 
soon be extracted by LHC experiments and exotic phenomena will indeed be observed. 
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